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ABSTRACT 

Ultraluminous supersoft X-ray sources (ULS) exhibit supersoft X-ray spectra with blackbody tem¬ 
peratures below 0.1 keV and bolometric luminosities above 10^^ ergs s“^. In this Letter we report the 
first optical spectroscopic observations of a ULS in M81 using the LRIS spectrograph on the Keck I 
telescope. The detected Balmer emission lines show the mean intrinsic velocity dispersion of 400 ± 

80 km s“^, which are consistent with from an accretion disk. The spectral index of the continuum 
on the blue side is also consistent with the multi-color disk model. The emission line exhibits a 
velocity of ^ 180 km/s relative to the local stellar environment, suggesting this ULS is possibly a halo 
system in M81 belonging to an old population. No significant shift is found for the Hq, emission line 
between two observations separated by four nights. 

Subject headings: galaxies: individual (M81) — X-rays: binaries 


1. INTRODUCTION 

Ultraluminous supersoft X-ray sources (ULSs) are 
pointlike, non-nuclear X-ray sources having extremely 
soft spectra with equivalent blackbody temperatures be¬ 
low 0.1 keV and bolometric luminosities above 10^^ 
ergs s“^. They are thought to be massive white dwarfs 
(WDs) burning accreted material on their surface or 
intermediate-mass black holes (IMBHs; 10^ — 10^ M<^ ) 
w ith sub-Eddington a ccretion (jLiu fc Di Stefanoll2QQ8^ . 

I Swartz et al.l (|2QQ2[ ) observed the nearby spiral galaxy 
M81 with Chandra ACIS, and discovered an intrigu¬ 
ing ULS in the bulge, R.A. = 09^55’^42.2^, dec. = 
+69°03'36.5" (J2000.0) (heareafter M81 ULSl). Its spec¬ 
trum can be fitted by a blackbody model with a temper¬ 
ature of ^ 70 eV. The bolometric luminosity calculated 
with the distance of M81 is ^ 10^^ ergs s“^. The follow¬ 
up X-ray study reveals that its spectrum can be described 
by either a blackbo dy for a W D or a multi-color accretion 
disk for an IMBH (|Liu Il200§ ). 

Its optical cou nterpart was detected b y Hubble Space 
Telescope (HST: [Liu fc Di Stefanol [200^ . The spectral 
energy distribution (SED) from broad-bands photome¬ 
try exhibits a blue and a red component. The spectral 
index of the blue component is consistent with a geomet¬ 
rically thin accretion disk, and the red component could 
be described by an AGB star. The SED also shows exces¬ 
sive Ha emission which is probably originated from the 
accretion disk or surrounding material photo-ionized by 
the soft X-ray emission of M81 ULSl, but the band pho¬ 
tometry results cannot provide the width of the emission 
line known as an indicator of physical process. 

The SED measurements, however, suffer from the in¬ 
trinsic variabilities of ULSl, since a flux decrease by a 
factor of ^ 2.4, as iTao et al l (j2Qllf ) presented, occurred 
in less than a week in optical wavelength. Spectral ob¬ 
servations, on the other hand, are independent of the 
variabilities. The expected Balmer emission lines in the 
spectrum will enable us not only to characterize its phys- 
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ical conditions but also to probe its local environments 
(|Moon et al.ll2QlIl) . This information is essential for us 
to understand the nature of M81 ULSl. 

In this Letter we report Keck spectroscopic observa¬ 
tions of M81 ULSl and present evidence for the existence 
of an accretion disk. Section [2] describes the Keck/LRIS 
observations and the result. The discussion is given in 
Section [3l 

2. DATA ANALYSIS AND RESULTS 

The observations of M81 ULSl occurred on 2010 April 
13 and 2010 April 17 during its expected X-ray-low state 
(jWang et al.l l2015f ) using the Low Resolution Imaging 
Spectrograph (LRIS) on the Keck 110m telescope. Three 
exposures of 1000 s were taken in the first night and two 
of 1200 s in the second night. The mean seeings in B 
band were 0.6" and 0.8" respectively. The light of the 
counterpart was masked with a 0.7" wide slit and split 
with a beam dichroic to the blue and red sides followed 
by a 300 lines mm“^ and a 400 lines mm“^ grating. 

The spectra were reduced in a standard way with 
IRAE. Eirst of all, raw EITS files were bias-subtracted, 
flat-corrected, and combined. On the blue side of the 
spectra, the position of the optical counterpart along the 
slit was verified by comparing its position in the target 
acquisi tion image (Eigure jH^) w ith an HST ACS E606W 
image (|Liu fc Di StefanQll2008f ). On the red side, since 
the counterpart was not obvious along the slit, an offset 
between the counterpart and the WD (PG1708 + 602) 
was used to verify the position of the counterpart. 

Subsequently, raw spectra of the counterpart were ex¬ 
tracted with an aperture size of 1". The wavelength cali¬ 
bration was then carried out based on the line lists given 
in the manual of Keclfl. The precision of the calibration 
is 0.2 A which is obtained from RMS (root mean square) 
of arc-lamp fitting. 

Einally, PG1708 + 602 was used as the standard star 
to calibrate the flux by app lying the standard flux given 
bv iMassev fc Strobell (|1990[ ). Since the standard flux of 

^ http://www2.keck.hawaii.edu/inst/lris/. 
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Fig. 1.— a: False-color Spitzer map of M81 used 5.8 p,m (blue)^ 
8.0 jum (green) and 24 j um (red). The blue line is 2" slit in 
IVega Beltran et aTI (|2001h . The light green rectangle presents the 
sky covered by Keck, b: The target acquisition image (gray) on 
the blue side and the image of 0.7" wide slit (red) in April 13. 
The green circle gives the position of the optical counterpart in 
a radius of 1". c: A section of the 2-D dispersed image in April 
13. The yellow arrows mark the locations of H^, Hq, and the un¬ 
known emission line around 5530 A. The red rectangle marks the 
Hq; absorption line in a region corresponding to the red rectangle 
in b. 

PG1708 + 602 covered the wavelength from 3126 A to 
8004 A, we had to extend the tabulated values to 10000 
A for the calibration on the red side (the upper panel in 
Figure [2]). 


2 . 1 . Baimer Emission Lines 
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Fig. 2.— Upper panel'. Observational result of PG1708+602. The 
fluxes are in the unit of 10“^^^ erg s“^ cm“^ A“^. Middle panel: 
The spectrum of optical counterpart. The integrated spectra of 
counterpart are smoothed to illustrate features clearly. Grey lines 
plot error spectra in 1 cr and the positions of hydrogen lines are 
given in red. The fluxes are in the unit of 10“^® erg s“^ cm“^ 
respectively. Lower panel: Emission lines of H^, Hq and the 
emission line around 5530 A after the subtraction of the continuum. 
Red lines stand for emission lines obtained in 2010 April 13, blue 
for lines obtained in 2010 April 17 and black for combined emission 
lines. The fluxes are in the unit of 10“^® erg s“^ cm“^ A“^. 


As shown in Figure [21 the Balmer emission lines are 
notable features in the spectra. Here we used min¬ 
imization to fit the centers and FWHMs (full width at 
half-maximum) of Ha and another notable emission 

line around 5530A. The results with Icr errors are listed 
in Table [H The mean FWHM of H /3 and Hq, emission 
lines derived from the fitting is 490 ± 80 km s“^, which 
is significantly larger than the spectral resolution, 280 
±10 km s“^measured from Hg I A5461 in the arc lamp 
spectrum. The intrinsic dispersion is 400 ± 80 km s“^for 
the Balmer emission lines. The radial velocity of ULSl 
is consistent with zero (an average of 2 ± 23 km s“^) 
within the precision of wavelength calibration, which im¬ 
plies that M81 ULSl is unlikely a distant AGN with a 
large receding velocity. 

The Balmer emission lines are mainly seen in spectra 
of H H regions, planetary nebulae and accretion disks 
around compact objects. The line dispersions for H H re¬ 
gions and planetary neb ulae should be comparable to in¬ 
strum ental dispersions ([Fang et ah! 120131: iNicholls et ah! 
120141) , while for accretion disks the line dispersions range 
from a few hundreds to thousands of kilo meters per sec¬ 
ond. The detected Balmer emission lines of M81 ULSl 
are significantly broader than the instrumental disper¬ 
sion, and are consistent with from an accretion disk 


around a compact object. 

The observations in two nights enable us to measure 
the shift of Hq related to the binary motion. We used 
the technology of the cross-correlation phasor to estimate 
the relative velocity. The technology is based on the 
correlation in the wave-number space and it can provide 
more information than the normal cross-corre lation, such 
as the significant level (see lMisra et al.ll 2 QlQl for detail). 
The shifts of Hq and sky emission lines are measured in 
the pixel space to avoid the uncertainty of the wavelength 
calibration. The relative shift of Hq to the sky line is 
0.43±0.19 pixel, corresponding to 0.61±0.38 A toward 
short wavelength. Then we used y-square minimization 
to check the relative shift and the result is 0.42±0.36 
pixel. No shift is found for Hq emission line at high 
significant level in the two nights observations with a 
span of 4 nights. 

2 . 2 . Balmer absorption lines 

Hq absorption features are detected in a specific region 
along the slit as illustrated in Figure [T]c. The regions with 
absorption features are located in the same area of the 
sky but different parts of the GGD in two nights observa¬ 
tions, so these absorption features are unlikely artifacts 
of the GGD. The same features were also found near H/ 3 , 
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TABLE 1 

Line fitting results 



4.13 

Center (A) 

4.17 

Combined 

4.13 

FWHM (A) 

4.17 

Combined 

~5530A 

m * 

^ cx,ah 

4861.4T1.1 

5532.3T1.5 

6562.9T0.2 

6558.8T0.5 

4861.2T0.3 

5543.1T1.6 

6562.8T0.2 

6558.9T0.3 

4861.3T0.6 

5536.6T1.6 

6562.8T0.2 

8.8T2.7 

32.5T3.5 

10.8T0.5 

8.2T0.7 

32.9T4.0 

9.4T0.3 

8.6T1.3 

32.2T3.8 

10.2T0.4 

Note. — 

* Calculated for the Ha 

absorption as marked by the red rectangle 

in Figlc. 




TABLE 2 

GALEX Observations and results for M81 ULSl 


Tile 

Band 

Effective 

Exptime 

Min 

Max 

Inner 

Outer 

AB Mag. 



Wavelength 


ObsDate 

ObsDate 

Radius 

Radius 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

GI1_071001_M81 

FUV 

1538.6 

14706.7 

2006-01-05 

2007-03-31 

6.3 

10.5 

22.27T0.10 


NUV 

2315.7 

29421.5 

2005-01-12 

2007-03-31 

7.9 

13.2 

21.36T0.17 


Note. — Col. (1): Tile name. Col. (2): Band name. Col. (3): Effective wavelength in angstrom. Col. (4): Total exposure time in 
seconds. Col. (5): Earliest observation U T date for visits which make up the coadd. Col. (6): Latest observation UT date. Col. (7): Inner 
radius in arcsecond for signal integration IMorrissev et HI (|2007I L Col. (8): Outer radius in arcsecond for background substruction. Col. 
(9): Magnitudes in AB system without correction of Calactic extinction. 



2000 4000 6000 800010000 
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Fig. 3.— SEP of M81 ULSl. T he blue and red components 
suggested bv ILiu &: Di Stefanol (120081 1 are plotted in red lines. Red 
points are photometri c results of GALEX and squares are those 
of HST obtained from ILiu fc Di Stefanol (120081 1. For correction of 
Galactic extinction, we h ave adopted the Gala ctic color excesses, 
E(B — E)=0.08, given bv ISchlegel et all (1199811 . an d the parame¬ 
teriza tion of the Galactic extinction law given by ICardelli et al.l 
(1198911 with an extinction ratio Ry = 3.1, and the conversion 
factors for the FUV and NUV are Apuv = 7.9E{B — V) and 
Ajsfuv — S.0E(B — V). 


but the signal is too low to derive reliable velocity. 

In order to investigate the association between the 
Balmer absorption lines and the environment of M81 
ULSl, false-color maps were constructed in optical 
{HST), UV {GALEX) and IR {Spitzer) wavelengths. 
Only the map of Spitzer exhibits obvious structures along 
the slit (Figure [T]a). The region with the absorption 
features is correlated with the dark area in the map 
of 8 jam, and the absorption features disappear when 


the 8 jam emission is strong. Note that there is a posi¬ 
tive correlation between 8 jam surface brigh tness and Hn. 
emiss ion powered by star forming activity (j Young et al.l 
l2Q14l1 . The observed absorption features appear in the 
non-star-forming regions, and likely come from the local 
stellar environments in the M81 bulge. 

The absorption features in the background exhibit high 
velocities, —180 ± 25 km s~^ at 52^^ from the nucleus . 
The stellar kinematics studies (jVega Beltran et ^l2QQl[1 
show that the bulge of M81 is rotating and the south-east 
side is approaching with velocities around —180 km s“^at 
45" from the nucleus. The similar velocity (—180 km s“^) 
suggests that the observed Hq, absorption lines come from 
stars in the bulge of M81. In comparison, the systemic 
velocity of M81 is —34 km s“^(obtained from NED). 

The emission line (—2 ± 23 km s exhibits a re¬ 
ceding velocity of 178 ± 35 km s“^relative to the ab¬ 
sorption line from the stellar environment. This suggests 
that ULSl is not a system comoving with stars in the 
bulge of M81. Such a relative velocity can come from an 
object in the halo of M81 projected onto its bulge, which 
is receding along the line of sight. Here we rule out the 
possibility that the system is located in the Milky Way, 
since the column density of neutral hydrogen atoms at 
the location of ULSl is 5.4 T 0 .5 x 10^^ cm“^ in the 
Milky Way (jGuver fc Ozell[2QQ^ . smaller than that de¬ 
rived from the X-ray spectrum fitting by a factor of 2 
^Liull2008ll . 

2.3. SED Construction 

The spectral observations of M81 ULSl enable us to 
obtain SED uncontaminated by its intrinsic variabili¬ 
ties. Here we removed all the emission lines and used 
a power-law function, /y oc A^, to fit the continuum of ^ 
4500—6000 A, the high response wavelength range, in or¬ 
der to estimate the spectral index on the blue side. The 
P derived from the best fitting is —2.36 ± 0.02, which 
corresponds to a ^ 0.36 ± 0.02 for oc The con- 
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tinuum on the blue side is consistent with the fj^ oc 
relation expected for the multi-color disk (MCD) model. 
On the red side, the spectrum of an AGB star modeled 
with a tempera ture of 2277 K and a rad ius of 1301 R©, 
as suggested bv lLiu fc Di Stefanol (|2QQ8[) . is presented in 
Figure [3l The combined SED is not consistent with the 
red-side spectrum, which suggests that the AGB model 
may be not a good explanation of the red component. 

As a connection of X-ray and optical wavebands, the 
UV emission is important to understand the SED of M81 
ULSl. In order to probe the UV emission of M81 ULSl, 
we used the arc hive data of the G alaxy Evolution Ex¬ 
plorer (GALEX [ iMdiitm et al.ll2QQ^ . M81 have been ob¬ 
served with three different tiles, and here we used the tile 
with exposure time over 10^ s to do accurate photometry. 
All the sub-exposures of the tile were taken during the 
expec ted low state of M81 ULSl (jLin II2QQ8I : IWang et all 
l2Q15f ) . The photometric results are listed in Table [2] and 
also shown in Eigure [3l Although with dispersions, the 
UV fluxes are lying along the MGD powerlaw. 

3. DISGUSSION 

In this letter we report the first optical spectroscopic 
confirmation of an accretion disk around a ULS. The 
broad Balmer emission lines of M81 ULSl revealed by 
the Keck/LRIS observations are consistent with from an 
accretion disk around a compact object. The spectral in¬ 
dex of the continuum on the blue side is also consistent 
with the fi; oc relation expected by MGD model. 
The velocity of the emission line relative to its local 
stellar environments suggests that M81 ULSl is possibly 
a halo system in M81 belonging to an old population. 
Gareful analysis shows no significant shift of Hq, emis¬ 
sion line with a span of 4 days. This might suggest a 
small inclination angle or a small velocity of the com¬ 
pact object. The latter can come from a long orbital pe¬ 
riod, or a very massive compact object such as an IMBH. 
However, the exact nature of the compact object is still 
unknown without monitoring of the long-term motion of 
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the system. 

Besides emission lines of Balmer series, an emission line 
arises around 5530 A on the blue side with an EWHM of 
1700 ± 200 km s“^, larger than that of Hq, by a factor of 
4 (the lower panel in Eigure [2]). All exposures obtained 
in two nights show the same emission feature and similar 
velocity dispersions, so it is unlikely an artifact or cos¬ 
mic rays. This emission line is too broad to be a nebular 
line, and it is probably not an Ee II emission line because 
Ee+ ion emits through a huge number of mnltiplets scat¬ 
tered across the blue side of th e spectrum ([Baldwin et al.l 
|2QQ4 iShapovalova et ani2Q12[ ) . Since the relative shift is 
10 A ± 2 A between two nights observations, the emis¬ 
sion line is unlikely related to the accretion disk. It is 
reported that N II A5530 and A5535 emission lines have 
been detected for Ae/Be stars (iMathew fc Bhavvall2QlQl : 

iMathew fc SubramanianJl2Qll[ ). However, no other fea¬ 

tures are found to support the existence of an Ae/Be star, 
such as Ee H or O I emission lines. More spectra with 
high signal-to-noise ratio are needed to draw a conclusive 
result. 

The spectrum of M81 ULSl shows a gap between the 
red and blue sides due to the low response on the edges 
of GGDs, and signal-to-noise ratio on the red side is very 
low. Eurther deeper spectroscopic observations with the 
coverage of ^ 5000—9000 A could present the complete 
spectrum with high signal-to-noise ratio on the red side, 
which will enable us to characterize the secondary. 
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